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ABSTRACT 

Microlensing is the only known direct method to measure the masses of stars 
that lack visible companions. In terms of microlensing observables, the mass 
is given by M = (c^/4G')fE^^E and so requires the measurement of both the 
angular Einstein radius, Be, and the projected Einstein radius, r^- Simultaneous 
measurement of these two parameters is extremely rare. Here we analyze OGLE- 
2003-BLG-238, a spectacularly bright (/mm = 10.3), high-magnification (Amax — 
170) microlensing event. Pronounced finite source effects permit a measurement 
of = 650 //as. Although the timescale of the event is only ^e = 38 days, one 
can still obtain weak constraints on the microlens parallax: 4.4 AU < te < 18 AU 
at the 1 a level. Together these two parameter measurements yield a range for 
the lens mass of 0.36 Mq < M < 1.48 Mq. As was the case for MACHO-LMC-5, 
the only other single star (apart from the Sun) whose mass has been determined 
from its gravitational effects, this estimate is rather crude. It does, however, 
demonstrate the viability of the technique. We also discuss future prospects for 
single-lens mass measurements. 

Subject headings: gravitational lensing — parallax 



1. Introduction 

When microlensing experiments were first proposed (Paczynski 1986) and implemented 
(Alcock et al. 1993; Aubourg et al. 1993; Udalski et al. 1993), it was not expected to be 
possible to measure the masses and distances of individual microlenses. The only microlens- 
ing parameter that depends on the mass and that is routinely measured is the Einstein 
timescale t-^, which is a degenerate combination of the lens mass M, and the lens-source 
relative parallax, TTrei, and proper motion, /irei- Specifically, 

o ^^TT 4G , mas , , 

iE = , ^E = V^MTTrel, « = ^— ~ 8.14— -, 1 

A*rel AU Mq 

where 9e is the angular Einstein radius. However, Gould (1992) showed that if both 9e and 
the microlens parallax. 
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could be measured, then the mass and lens-source relative parallax could both be determined, 

M=^, 7rrel = 7rE^E. (3) 

Nevertheless, of the roughly 2000 microlensing events detected to date, there have been 
only of order a dozen for which 6'e has been measured and a dozen for which tte has been 
measured. Moreover, there is only one event, EROS-BLG-2000-5, with measurements of both 
parameters, and so for which the microlens mass and distance have been reliably determined 
(An et al. 2002). Since this one event was a binary, and since all the other stars with directly 
measured masses are components of binaries, it remains the case today that the only single 
star with a directly measured mass is the Sun. 

The one partial exception is the microlens in MACHO-LMC-5. Alcock et al. (2001) were 
able to measure both 9-e and tte for this event and so measure the mass and distance. These 
estimates were completely inconsistent with photometry-based estimates of these quantities, 
but Gould (2004) resolved this puzzle by showing that the tte measurement was subject to a 
discrete degeneracy for this event and that the alternate solution was consistent at the few 
a level with the photometric evidence. Nevertheless, since the error in the mass estimate is 
about 35%, this mass determination must be regarded as very approximate. 

Here wc analyze OGLE-2003-BLG-238, the brightest microlensing event ever observed 
and only the fourth reported point-lens (i.e., non-binary) event with pronounced finite-source 
effects. As with the other three such events (Alcock et al. 1997; Smith, Mao & Wozniak 
2003; Yoo et al. 2004), these finite-source effects allow one to measure ^e with reasonably 
good (~ 10%) precision, where the error is typically dominated by the modeling of the source 
rather than the microlensing event. Hence, if tte could also be measured, it would be possible 
to determine M. 

Despite the event's short duration, it is still possible to detect parallax effects in OGLE- 
2003-BLG-238 because of its bright source and high magnification. For short events hke 
this one, the Earth's acceleration can be approximated as uniform during the event. Gould, 
Miralda-Escude & Bahcall (1994) showed that under these conditions, the parallax effect 
reduces to a simple asymmetry in the lightcurve around the peak. The high magnification of 
OGLE-2003-BLG-238 permits a very accurate measurement of the peak time of the event, 
which in turn makes the fitting process very sensitive to this small asymmetry. The brightness 
of the source allows high precision photometric measurements even in the wings of the 
event, which enable detection of these subtle deviations. Unfortunately, as also shown by 
Gould et al. (1994), the simphcity of the parallax effect for short events imphes that only 1- 
dimensional parallax information can be effectively extracted, whereas the microlens parallax 
is intrinsically a 2-dimensional vector tte. That is, while one component of the vector parallax 
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is well determined, the scalar parallax tte is not well determined, and this degrades the mass 
determination through equation (3). Nevertheless, this is still only the second single star 
(other than the Sun) for which any direct mass measurement at all can be made. 

2. Observational Data 

The microlensing event OGLE-2003-BLG-238 was identified by the OGLE-III Early 
Warning System (EWS) (Udalski 2003) on 2003 June 22. It peaked on HJD' = HJD-2450000 = 
2878.38 (Aug 26.88) over South Africa. OGLE-III observations were carried out in / band 
using the 1.3-m Warsaw telescope at the Las Campanas Observatory, Chile, which is op- 
erated by the Carnegie Institution of Washington. While OGLE-III normally operates in 
survey mode, cycling through the observed fields typically once per two nights during the 
main part of the bulge season, it can switch rapidly to follow-up mode if an event is of partic- 
ular interest and requires dense sampling. The high magnification of OGLE-2003-BLG-238, 
which was predicted while the event was developing, and possible deviations from a single- 
point-mass microlensing lightcurve profile were the main reasons that OGLE observed this 
event in follow-up mode. 

The OGLE-III data comprise a total of 205 observations in I band, including 144 during 
the 2003 season and 61 in the two previous seasons, 2001 and 2002. The exposures were 
generally the standard 120 s, except for three nights (HJD' 2877.5-2879.7) around the maxi- 
mum when the star was too bright for the standard exposure time. The time of exposure was 
adjusted then to the current brightness of the lens and seeing conditions to avoid saturation 
of images and was as short as 10 s on the night of maximum. Photometry was obtained 
with the OGLE-III image subtraction technique data pipeline Udalski (2003) based in part 
on the Wozniak (2000) DIA implementation. 

Following the alert, the event was monitored by the Microlensing Follow-Up Network 
(yuFUN, Yoo et al. 2004) from sites in Chile and Israel, and by the Probing Lcnsing Anomalies 
Network (PLANET, Albrow ct al. 1998) from sites in Chile and Tasmania. The /iFUN Chile 
observations were carried out at the 1.3m (ex-2MASS) telescope at Cerro Tololo InterAmer- 
ican Observatory, using the ANDICAM camera, which simultaneously images at optical and 
infrared wavelengths (DePoy et al. 2003). There were a total 203 images in I from HJD' 
2870.5 (eight days before peak) until HJD' 2950.5 at the end of the season. The exposures 
were generally 300 s, but were shortened to 120 s for 81 exposures over the peak. The 
exposures should have been further shortened on the night of the peak but, due to human 
error, this did not happen. Hence, all 19 of these images were saturated. yuFUN obtained 12 
points in V , primarily to determine the source color. This includes one saturated point over 
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the peak. All saturated images were discarded. 

The //FUN Israel observations were carried out in / band using the Wise Im telescope at 
Mitzpe Ramon, 200 km south of Tel- Aviv. There were 14 observations in total, all restricted 
to the peak of the event, 2877.3 < HJD' < 2883.3. The exposures (all 240 s) were obtained 
using the Wise Tektronix IK CCD camera. All //FUN photometry was extracted using 
DoPhot (Schechter, Mateo & Saha 1993). 

The PLANET Chile observations were carried out in R band using the Danish 1.54m 
telescope at the European Southern Observatory in La Silla, Chile. There were a total of 
68 observations from HJD' 2874.6 to HJD' 2883.7. The exposure times ranged from 2 to 80 
seconds. The PLANET Tasmania observations were carried out in / band using the Canopus 
Im telescope near Hobart, Tasmania with 52 observations from HJD' 2877.9 to HJD' 2903.9. 
The exposure times ranged from 60 to 300 seconds. During the first night in Tasmania the 
data were taken despite significant cloud cover by observing through "holes" in the cloudy 
sky. This proved feasible only because of the extreme brightness (/ ~ 11) of the source 
and demonstrates the importance of carefully monitoring events in real time to determine 
whether they should be observed despite truly awful conditions. 

The position of the source is R.A. = 17M5°^50".34, decl. = -22°40'58'.'l (J2000) (/, h = 
5.72, 2.60), and so was accessible for most of the night near peak from Chile and Tasmania, 
but for only a few hours from Israel. The combined data sets are shown in Figure 1. 



In general, the fiuxes, Fi{t), observed during a microlensing event by i = 1, . . . , n obser- 
vatory/filter combinations are fit to the form. 



where j is the fiux of the unmagnified source star as seen by the ith observatory and F^, j 
is any background flux that lies in the aperture but is not participating in the microlensing 
event. (The one exception would be a binary-source event, in which case two source-star 
terms and two magnification functions would be required.) 

In most cases, the lensed star can be fairly approximated as a point source. The mag- 
nification is then given by (Paczynski 1986), 



3. 



Finite-Source Effects 



(4) 



A{u) 



(5) 



li(li2 + 4)1/2' 
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where u is the angular source- lens separation in units of the angular Einstein radius ^e- 
However, this approximation breaks down for u < p, where, 

P - ^ (6) 

is the angular source size normalized to 6e- Finite-source effects then dominate. An 
appropriate formalism for incorporating these effects is given by Yoo et al. (2004). Here we 
summarize the essentials. If limb darkening is neglected, the total magnification becomes 
(Gould 1994; Witt & Mao 1994; Yoo et al. 2004), 

A,^,{u\p) ~ A{u)Bo{u/p), Bo{z) = -zE{k, z), (7) 

TT 

where E is the elliptic integral of the second kind and k = min(2;~^,l). This formula is 
accurate to 0(p^/8) (Yoo et al. 2004), and hence it applies whenever (as in the present case) 
p <^1. Note from Figure 1 that the finite-source fit passes first above the point-source curve 
and then moves below it. This transition occurs when Bo{z) — 1, which (from fig. 3 of Yoo et 
al. 2004) occurs at z ~ 0.54. By contrast, for MACHO-95-30 (Alcock et al. 1997) and OGLE- 
2003-BLG-262 (Yoo et al. 2004), the finite-source fits remain above the point-source curves 
throughout because in those cases the minimum source-lens separation (impact parameters) 
were Uq ~ 0.7p and Uq ~ 0.6p, respectively. 

To include the effects of hmb darkening, we model the source profile Sx at each wave- 
length A by, 

Sx{'d)^Sx[l-rx{l-^cos^)], (8) 

where Tx is the linear limb-darkening coefficient and is the angle between the normal to 
the stellar surface and the line of sight. The magnification is then given by, 

Aid(«|p) = A{u)[Bo{z) - rB,{z)], (9) 

where -61(2;) is a function described by equation (16) and figure 3 of Yoo et al. (2004). 



4. Parallax Effect 

If the motions of the source, lens, and observer can all be approximated as rectilinear, 
the source-lens separation, u, in equation (5), can be written, 

<t) = v[r{tw + mm (10) 
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where, 

r(t) = ^, m-uo. (11) 

The simplest point-source fit to a microlcnsing event requires five parameters, the source 
flux Fg, the background flux F^, the time of closest approach to, the Einstein time scale t^, 
and impact parameter uq. 

However, even if the source and lens can be assumed to be in rectilinear motion, the 
Earth is not. Especially for the long events (ts > yr/27r), the Earth parallax effect must be 
taken into account. The event OGLE-2003-BLG-238 lasted only 38 days, and hence parallax 
effects would be negligible if the source were not very bright and highly magnifled, both of 
which facihtate detection of the very subtle parallax effect. Moreover, after it was realized 
that flnite-source effects had been detected, both OGLE and //FUN intensifled observations 
of the event in the hope of measuring the parallax and so combining the result with the 
source-size measurement to determine a mass. 

Historically, parallax fits were carried out in the heliocentric frame. That is, uo was 
adopted as the point of closest approach to the Sun, and to was the time at which this 
approach occurred. When, as in the present case, parallax is only weakly detected, the 
trajectory relative to the Sun is poorly determined, so to and Uq have very large errors that 
are highly correlated with the parallax parameters. In the geocentric frame, by contrast, 
to and Mo are directly determined from the time and height of the peak of the light curve 
(Dominik 1998). Gould (2004) further refined the geocentric frame by subtracting out the 
Earth-Sun relative velocity as well as their positional offset. This frame is appropriate for the 
analysis of OGLE-2003-BLG-238. Hence, we follow the Gould (2004) formalism for modeling 
parallax in the geocentric frame. 

The parallax effect is parameterized by a vector tte whose magnitude gives the ratio of 
the Earth's orbit (1 AU) to the size of the Einstein ring projected onto the observer plane 
(ffi) and whose direction is that of the lens-source relative motion as seen from the Earth at 
the peak of the event. Exphcitly, tte = AU/fE. Equation (11) is then replaced by, 

r{t)^'^-^ + ST, m= ^0 + 5(3, (12) 

IE 

where, 

{St, 5(5) = TTE As = (tte • As, tte x As), (13) 

and As is the apparent position of the Sun relative to what it would be if the Earth remained 
in rectilinear motion with the velocity it had at the peak of the event. See equations (4)- (6) 
of Gould (2004). More exphcitly, 

{5t, 5(3) = [AsN(t)7rE,N + AsB(i)7rE,E, -AsN(t)7rE,E + AsE(t)7rE,N], (14) 
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where the subscripts N and E refer to components projected on the sky in north and east 
celestial coordinates. 

A major advantage of this formalism is that the parameters to, Uq, and Ie are virtually 
the same for the parallax and non-parallax solutions and, especially important in the present 
case, when the parallax solution is varied in the tte plane to obtain likelihood contours. 

4.1. NonpciralleLx Fit 

As we will show in § 4.2, the parallax effect is quite subtle and so can, to a first 
approximation, be ignored. On the other hand, the finite-source effects are quite severe 
(see Fig. 1). We therefore begin by fitting the hghtcurve by taking into account finite- 
source effects (including linear limb darkening) but not parallax. The fit therefore contains 

6 geometric parameters (tg, "^^oi ^e? T/, Tr) as well as 12 flux parameters (two for each 
observatory/filter combination). The results are listed in Table 1 and plotted in Figure 1. 
Also shown in Figure 1 is the lightcurve that would have been generated by the same event 
but assuming that the source had been a point of light. As discussed in § 3 this remains 
below the finite-source curve until z = u/p = 0.54, and then rises dramatically above it. 

4.2. Parallax Fit 

To find the best-fit parallax tte and the error ellipse around it, we conduct a grid 
search over the tte plane. That is, we minimize by holding each trial parameter pair 
tte = (tte.at, tte.e) fixed, while allowing the remaining 18 parameters (see § 4.1) to vary. 
After finding the best fit tte, we hold this fixed and renormalize the errors so that P^r 
degree of freedom is unity. We eliminate the largest outlier point and repeat the process 
until there are no 3 cr outhers. Of course, we first ehminate the 19 saturated points in //FUN 
Chile / and 1 saturated point in //FUN V. We find that this procedure removes 5 points 
from the OGLE data, an additional 17 points from the yuFUN Chile I data, 7 points from 
the PLANET R data, 2 points from the PLANET / data, and none from the other data 
sets. The final renormalization factors are 1.96, 0.87, 2.2 and 1.4 for OGLE, yuFUN Chile 
/, PLANET Tasmania / and PLANET Chile R data, respectively. The other two data 
sets, //FUN Chile V and /iFUN Israel /, do not require renormalization. This cleaned and 
renormalization data set is used in all fits reported in this paper and is shown in Figure 1. 
It contains 200 points from OGLE I, 167 from /tFUN Chile I, 14 from /(FUN Israel 7, 11 
from /tFUN Chile V, 50 from PLANET Tasmania I, and 61 from PLANET Chile R. 
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Figure 2 shows the resulting hkehhood contours in the tte plane. The best fit is at 
(tte.e, 7rE,Ar) = (0.0664,-0.0205). The contours are extremely elongated with their major 
axes almost perfectly aligned with the North-South axis. Gould et al. (1994) showed that 
short events would yield essentially 1-dimensional parallax information because the Earth's 
acceleration vector is basically constant over the duration of the event. Hence, only a single 
parallax parameter can be measured robustly, namely the magnitude of the asymmetry of 
the lightcurve. This yields information about the component of the projected lens-source 
relative velocity parallel to the Earth's acceleration (projected onto the plane of the sky) but 
not the perpendicular component. At the peak of the event, the position of the Sun (pro- 
jected onto the plane of the sky - see eq. [6] of Gould 2004) is {se, sn) = (-0.930, 0.028)AU, 
which means that the Earth is accelerating in the same direction. Hence, one expects the 
direction of maximum sensitivity (minor axis of the error ellipse) to be at a postition angle 
tan^-'^(— 0.930/0.028) = 91°. 725 (north through east), which agrees quite well with the ori- 
entation (91°.769) shown in Figure 2. The event MOA-2003-BLG-37, which was also short 
(^E ~ 42 days) shows similar highly elongated parallax-error contours (Park et al. 2004). 

Figures 1 and 2 both illustrate that the parallax effect in OGLE-2003-BLG-238 is weak. 
The residuals in Figure 1, which shows the fit without parallax, demonstrate that the asym- 
metry is quite subtle. Figure 2 shows that the error contours extend almost to the origin. 
That is, from Table 1, the addition of two parallax parameters reduces from 510.6 to 
478.3, a 5.5 a effect. 

4.3. Check for Pcirallctx Degeneracies 

Gould (2004) showed that microlensing events, particularly those with short timescales 
(^E ^ yi'/27r), could be subject to a discrete four-fold degeneracy. One pair of degener- 
ate solutions, which was previously discovered by Smith, Mao Sz Paczyhski (2003), takes 
Uq —Uq. The remaining parameters are then similar to the original parameters, with 
the differences being proportional to Uq. Since in the present case Uq is extremely small, 
uo — 2 X 10~^, one expects that these two solutions would be virtually identical, and this 
proves to be the case. 

The other pair of solutions arises from the jerk-parallax degeneracy, which predicts that 
if tte = (tte,!!, tte.x) is a solution, then tt^ = (tt^ y, tt^ j^) is also a solution, with 



^E,|| = tte,!!, tTe,^ = -(7rE,± + tt^-x). 



(15) 
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where TVj is the "jerk parallax" . In the approximation that the Earth's orbit is circular, 

^_4j^r sin/3ec , . 

""''^ 3 27rtE (cos2 sin^ + sin^ V)^/' ' 

where /9cc is the ecliptic latitude of the source and ip = 69° is the phase of the Earth's 
orbit relative to opposition at the peak of the event. In the present case, the event is quite 
close to the ecliptic, Pec ~ 0°.8, so tTj^j, ~ —0.037. Since tte,± = 0.018, this implies that 
TTg = ~(^E,± + '^j,±) ~ 0.017, which is almost identical to tte.x- Hence, no degeneracy is 
predicted, and this expectation is confirmed by Figure 2, which shows a single minimum. 

Note that for events seen right on the ecliptic, /?cc = 0, equation (15) predicts 71"^^ = 
— 7rE,±. Indeed, for this special case, the degeneracy is exact to all orders and not only to 
fourth order in the perturbative expansion as was derived by Gould (2004). That is, since the 
accelerated motion is exactly along the ecliptic, there is no way to distinguish whether the 
component of lens-source relative motion perpendicular to the ecliptic is toward the north 
or south. Since OGLE-2003-BLG-238 is very near the ecliptic, one would naively expect it 
to be strongly affected by this degeneracy. In fact, it is only because 7rE,± is also very close 
to zero that the degeneracy is avoided. 

The extreme axis ratio of the parallax-error ellipse, c(7rE,±)/(j(7rE,||) = 17, confirms that 
the information about parallax is essentially 1-dimensional, as predicted by Gould et al. 
(1994). However, it is not perfectly 1-dimensional: while tte.j. is highly consistent with zero 
there is some constraint, however weak, on this quantity. We search for the origin of this 
constraint within the context of the Gould (2004) formalism. For the limiting case (relevant 
here) of -Uq ~^ 0, 7rE,± first enters in the fourth order term C4 in the Taylor expansion 
— ^iCi{t — toY- For sufficiently large 7rE,_L, C4 ~ (Q;7rE,±/2)^, where a ~ (60days)~^ 
is the apparent acceleration of the Sun at the peak of the event, divided by an AU (see 
eq. [20] of Gould 2004). Hence, for uq (i.e., essentially everywhere in the present case), 
C^lt — to)^ ~ 0.04M*^7r| j^, implying a perturbation Au ~ 0.02ii^7r| j^. If we now consider 
tte.x = 0.4 (roughly the 2a upper limit), and focus on -u ~ 1.5 (where there is still a high 
density of 1% photometry and where dlnA/du ~ 0.3 is still fairly high), then the amplitude 
of the effect is small, AA/A ~ d\nA/du x 0.02M^7r|_|_ ~ 3 x 10"^, but stiU plausibly large 
enough given the large number of relatively high precision measurements. 

5. Negative Blending 

As seen from Table 1, the best nonparallax fit of the background flux for OGLE / 
band, /xFUN Chile / band and V band, and PLANET / band are all negative. There are 
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three potential reasons for such negative background fluxes: systematic photometry errors, 
unmodeled effects in the hghtcurve that are absorbed by the blended flux parameter, and 
"negative flux" from unlensed sources. The flrst possibility is virtually ruled out by the fact 
that the negative background flux appears in so many unassociated lightcurves. The last 
possibility is not as ridiculous as it might flrst appear because the dense Galactic bulge flelds 
have a mottled background of turnoff and main-sequence stars. If the source happens to lie in 
a hole in this background, it will appear as negative (Park et al. 2004). The -Fb/F^ = —5% 
from the OGLE photometry (which has the most extensive basehne) , would require a "hole" 
corresponding to a star 20 times fainter than the source i.e, /o,"hoie" ~ 17.6 (see Fig. 3). This 
is a plausible brightness for a hole in the unresolved turnoff stars. Combining this value 
with the Fb/Fs = —17% measurement from PLANET R, yields a color difference between 
the "hole" and the source, 

A(R -I) = {R- -{R- /), = -1.25 ± 0.20, (17) 

whereas the expected value (given the source position in Fig. 3) is about A{R — /) ~ —0.5. 
Hence this explanation is not completely self-consistent. 

Because the effect of the blending parameter is even about the peak, it can absorb effects 
of other even parameters including Fg, p, uq, and ^e- Since all of these are taken into account 
in the nonparallax flt (and its errors) these cannot be the cause. However, as pointed out 
by Smith et al. (2003) microlcns parallax can also mimic blending. Within the formalism of 
Gould (2004), the blending fraction is correlated with 7Te,± which is also an even parameter 
(see Park et al. 2004). 

The best-flt parallax solution still contains negative background fluxes, although these 
are slightly reduced in magnitude relative to the nonparallax flt, while the errors are some- 
what increased. The reduction reflects the absorption of some of the negative blending into 
7rE,±, while the larger errors reflect the covariance between Fj, and 71e,±- However, since the 
negative blending is still detected with substantial significance, parallax cannot be the whole 
story. A "hole" in the mottled bulge background of turnoff stars remains the most plausible 
explanation for the negative blending, although as discussed above, this explanation is not 
perfect. 



6. Error Determination 

We use Newton's method to find the minimum with respect to the 18 parameters of 
the nonparallax model. This procedure utilizes the Fischer matrix, and therefore automati- 
cally generates a covariance matrix and so linearized error estimates for all the parameters. 
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We find, however, that Newton's method fails when we add the two parallax parameters, 
probably because the effect is too subtle to withstand the numerical noise induced by nu- 
merical differentiation of the finite source effects. We therefore hold tte fixed at a grid of 
values and, at each one, minimize with respect to the remaining 18 parameters. The re- 
sulting contours are shown in Figure 2. To estimate the errors we use the method of "hybrid 
statistical errors" given in Appendix D of An et al. (2002). First, Newton's method auto- 
matically yields Qj, the covariance matrix of the model parameters (collectively Oj) with the 
two parameters tte held fixed at their best-fit values. Next we evaluate the two-dimensional 
covariance matrix c^„, where m, n range over the parameters (tte.at, tte.b), by fitting the 
contours in Figure 2 to a parabola. Third we evaluate dai/dam, the change in the best- fit 
model parameter as one of the parallax parameters is varied over the grid. Finally, we 
evaluate the covariance matrix by. 



We determine the angular size of the source from the instrumental color-magnitude 
diagram (CMD), using the method developed by Albrow et al. (2000b) and references therein, 
which is concisely summarized by Yoo et al. (2004). We measure the offsets in color and 
magnitude between the unmagnified source star and the center of the clump giants, A/ = 
Is - /clump = -0.02, A(\/ -I) = {V - I)s -{V - /)ciump = 0.22. For the dereddcncd color 
and magnitude of the clump center, wc adopt [{V — /)o, /o]ciump = (1.00, 14.32) (Yoo et al. 
2004), then transform from (V" — /)o to (V — K)q using the color-color relation of Bessell 
& Brett (1988). We obtain [{V - /)o,/o]. = (1-22, 14.30). Finally, using the color/surf ace- 
brightness relation of van Belle (1999), we obtain = 8.35 ± 0.72 //as, where the error is 
dominated by the 8.7% scatter in the van Belle (1999) relation. 

From the best-fit value p — 0.0128, we then obtain. 



= 652 ± 56 //as, /(,ei = 6.20 ± 0.54 mas yr"^ = 29.4 ± 2.6 km s"^ kpc~\ (19) 




(18) 



7. 



Estimates of Meiss and Distance 



7.1. Measurement oi 6e 
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7.2. Mass and Distance Estimates 

The best parallax fit for the event is tte = 0.0695, which when combined with equa- 
tions (3) and (19) yields, 

M = 1.15Mq (best fit). (20) 

However, the errors are quite large. Figure 2 shows that at the 1 a level, tte lies in the range 
0.2256 > tte > 0.0552, which implies 

0.36 < M/Mq < 1.48 (1 a). (21) 

The same microlens parallax estimates lead (through eq. [3]) to a best relative parallax 
estimate of TTrei = 45 //as and a 1 cr range of 147 //as > TTrei > 35 /tas. If one adopts a source 
distance of — 8kpc, this corresponds to a distance range 3.6 kpc < A < 6.3 kpc. 

At the 2a level, 0.4180 > tte > 0.0434, which leads to a mass range 0.19 < M/Mq < 1.86 
and a relative parallax range 273 /tas > TTrei > 28 /fas, corresponding to 2.5 kpc < Di < 
6.5 kpc. Therefore, the 2 a interval is consistent with most of the stellar range, but it does 
not provide any "new" information about the lens other than ruling out stellar-mass black 
holes and very late-type M dwarfs and brown dwarfs. On the other hand, it does serve 
as a basic consistency check on the viability of the microlens mass measurements, since if 
the method were plagued by strong systematic errors one would not necessarily expect the 
estimated mass to lie in the stellar range. 



7.3. A Single Star? 

As noted in § 1, part of the motivation for microlensing mass measurements is that 
this is the only known way to directly measure the mass of single stars. But how confident 
can we be that OGLE-2003-BLG-238 is a single star? In a future paper, we will analyze 
the lightcurve for the presence of all types of companions to the lens, particularly planetary 
companions. However, by simple inspection of the lightcurve, it is possible to place rough 
limits on stellar-mass companions, i.e., those with mass ratios q > 0.1. No such companions 
are possible with separations (in units of 9-e of the observed lens) of d ~ 1, since there would 
be clearly visible caustic crossings. 

As the separation is increased, the magnification pattern approaches a Chang-Refsdal 
lens with shear 7 = q.w/d'^ (Chang & Refsdal 1979, 1984). The width of the Chang-Refsdal 
caustic in the limit of 7 <g 1 is £ ~ 47. Prom inspection of the hghtcurve and the fact that 
1*0 ~ p/6 (see Table 1), we can say that such a caustic would certainly have been noticed if 
£ > p/2. This implies a hmit > {^qw/py^ = 25^^^^ 
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According to the so-called d — > duality discovered by Dominik (1999) and further 
elaborated by Albrow et al. (2002), the central caustics of extremely close binaries mimic 
those of wide binaries. Keeping to the convention that the observed Einstein radius cor- 
reponds to the mass near the observed peak in the lightcurve (i.e., the combined mass of 
the binary in the close case but just the mass one star in the wide case), one finds that 
7 = dlqc/{l + ^c)^- Thus, by the same argument as above, dc < 0.04(1 + Qcjqc^^'^- Combin- 
ing these two arguments and making use of equation (19), implies that the entire range of 
companion projected separations r±, 

0.2 AU (g^^ + g-^/2) ^ < < 85 AU qH^ ^ (excluded companions) , (22) 

is excluded. Here Di is the distance to the lens and Rq — 8kpc. Hence, if the lens has a 
stellar companion, it is either extremely close or very far away. 



8. Future Prospects 

To date, microlensing mass measurements of single stars have depended on very rare 
combinations of circumstances. The problem is somewhat more severe than was outlined in 
§ 1. There we noted that tte and 9e were separately measured only infrequently, so combined 
measurements are even more infrequent. However, for single stars, the events most likely 
to show the parallax effects from which one could measure tte are also the least hkely to 
show the finite-source effects from which one could measure ^e- That is, microlens parallax 
measurements generally require long events, > yr/27r, which tends to favor large masses, 
since tE oc M^/^. But the probability of significant finite-source effects (i.e., uq < p) scales 
as p = 9.^/9e, which favors small masses, since 9e oc M^/^. The combination of large ^e and 
small 9e implies low relative proper motion fi^-d = ^^e/^e- 

Actually, neither of the two single-star events with microlens mass measurements meets 
this criterion. OGLE-2003-BLG-238 has fi ~ 6masyr~^, which is a typical value for disk 
lenses seen against the bulge and is substantially higher than the proper motions expected 
for bulge-bulge events. MACHO-LMC-5 has /i ~ 20masyr~^, which of course is extremely 
fast. What can be learned about the future prospects for single-star mass measurements 
from the failure of both of these events to "fit the mold" ? 

OGLE-2003-BLG-238 was not long enough to obtain a good microlens parallax mea- 
surement. That is, only the tteji component of the microlens parallax vector tte can really 
be said to have been "measured". The other (7rE,±) component was only grossly constrained. 
Regarding MACHO-LMC-5, it was neither long enough for a very accurate measurement 
of tte, nor did it exhibit the finite source effects that are normally required to measure ^e- 
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The fact is that auxiliary, non-microlensing, data were needed to measure for this event. 
That is, the source and lens were separately resolved six years after the event, and from the 
measurement of their separation, Alcock et al. (2001) were able to deduce /Xj-ei, which when 
combined with microlensing data yielded 9e- 

This experience with MACHO-LMC-5 points to one future route to microlens mass 
measurements: give up altogether on measuring 9^ from the microlensing events: just focus 
on long events with good parallax measurements and measure 6-^ from post-event astromctry. 
Han & Chang (2003) estimated that 22% of disk-bulge events could be resolved 10 years after 
the event, assuming a resolution of 50 mas. 

Alcock et al. (2001) proposed a second route to measuring the mass of MACHO-LMC-5: 
use the fact that the source was already resolved to measure both the lens-source relative 
parallax tTj-ci and the lens-source relative proper motion fij-^. One sees directly from equa- 
tion (1) that if these two measurements are combined with a measurement of ^e, one obtains 
the lens mass M. Indeed, this would be a variant of the original idea of Refsdal (1964) to 
determine single-star masses by first finding nearby stars passing close to the line of sight of 
more distant stars and by then obtaining TTrei, and the angular deflection A^rei the impact 
parameter b — uo9e, all from astrometry. The only difference being that for these nearby 
lenses, which generally pass well outside the Einstein ring (uq ^ 1), 9e = {bA9rciY^'^ is 
determined directly from astromctry, rather than from the combination of astrometric (//rei) 
and photometric (^e) parameters employed by Alcock et al. (2001). 

Yet a third route is suggested by the experience of OGLE-2003-BLG-238. In spite of 
its short duration, the microlens parallax is well measured in one direction. If the lens 
could be resolved by post-event imaging, this would not only yield the magnitude of the 
proper motion fij.^i, but also its direction. Since the directions of and tte are the same, 
the proper motion measurement would at the same time resolve the parallax degeneracy. 
It may prove difficult to apply this method to OGLE-2003-BLG-238 itself because (from 
Table 1), the source is so much brighter than the lens. However, Ghosh et al. (2004) argue 
that another event, OGLE-2003-BLG-175/MOA-2003-BLG-45, shows excellent promise for 
yielding a mass with this method. 

Finally, a fourth route has been proposed by Agol et al. (2002). For very massive lenses, 
primarily black holes, the events will generally be long enough to measure tte, while 9-e may 
be large enough to measure it using precise astrometry from the deviation of the centroid 
of lensed hght relative to the source position. F. Delplancke (2004, private communication) 
expects that she and her team at the Very Large Telescope will achieve the required high 
precision for K < 11, 13, and 16 sources in 2004, 2005, and 2006, respectively. 
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In the longer term, it should be possible to measure single-star masses using the Space 
Interferometry Mission (SIM) in two distinct ways. First, SIMs high (4 //as) precision will 
allow one to carry out the original Refsdal (1964) proposal, provided appropriate lens-source 
pairs can be found (Paczynski 1998; Gould 2000; Salim & Gould 2000). Second, for mi- 
crolensing events generated by distant lenses (whether luminous or not) and for sufficiently 
bright sources, it will be possible for SIM to routinely measure 9e from the centroid dis- 
placement discussed above (Boden, Shao & Van Buren 1998; Paczynski 1998). On the other 
hand, since SIM will be in solar orbit, comparison of SIM and ground-based photometry 
of the event will yield microlens parallaxes according to the method of Refsdal (1966) and 
Gould (1995). Of order 200 such measurements should be feasible with the 1500 hours of 
5'/Mtime that has been allocated to this project (Gould & Sahm 1999). 
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Table 1. OGLE-2003-BLG-238 Fit Parameters 



Parallax Fit Nonparallax Fit 



Parameters Value Error Vaule Error 



to (days) 


2878.38123 


0.00078 


2878.38026 


0, 


.00076 


Uo 


0.00200 


0.00021 


0.00222 


0, 


.00019 


tE(days) 


38.18743 


0.22142 


37.58892 


0, 


.18946 


P 


0.01282 


0.00012 


0.01299 


0, 


.00012 


P/ 


0.47696 


0.06007 


0.46658 


0, 


,05793 


Pfl 


0.53287 


0.08419 


0.52438 


0, 


.08308 


tte.n 


- 0.02053 


0.19697 


0.0 








0.06639 


0.01328 


0.0 








0.06700 


0.01181 


0.0 








0.01847 


0.19706 


0.0 








-0.03716 


0.00763 


-0.05265 


0, 


.00549 


iFb/Fs)i, 


0.00959 


0.00704 


-0.01235 


0, 


.00479 


{Fb/Fs)v, 


- 0.00594 


0.10906 


-0.02527 


0, 


.07587 


iFb/Fs)i, 


0.02762 


0.08100 


0.02850 


0, 


.05617 


{FblFs)u 


-0.07444 


0.02112 


-0.09285 


0, 


.01485 


{Fb/Fs)R, 

x' 


-0.14948 


0.04248 


-0.16968 


0, 


.02990 


478.323 




510.643 







Note. — Observatories: 1=0GLE, 2=;uFUN Chile, 3=mFUN 
Israel, 4=PLANET; Ax^ :^ 32 




2870 2880 2890 

HJD-2,450,000 



Fig. 1. — Photometry of microlensing event OGLE-2003-BLG-238 near its peak on 2003 
August 26.88 (HJD 2452878.38). Data points are in / (OGLE: black; //FUN Chile: red; 
PLANET Tasmania: magenta; /iFUN Israel: green), V (/xFUN: cyan) and R (PLANET 
Chile: yellow). The saturated and other excluded points are not shown. The circles are 
displayed at different sizes to make the figure more readable. All bands are linearly rescaled 
so that Fs and are the same as the OGLE observations, which define the magnitude scale. 
The solid curve shows the best-fit nonparallax model with finite-source effect. The dashed 
curve shows the lightcurve expected for the same lens model, but with a point source. 
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Fig. 2. — Likelihood contours (cr = 1, 2, 3, 4, 5) as a function of the vector parallax tte- 7rE,|| 
gives the direction of the Sun's apparent acceleration. As expected from theory, the parallax 
is well constrained in this direction but poorly constrained in the orthogonal (7rE,±) direction, 
which hes at an angle of 1°.769 from the North-South axis (vertical line segment). 
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Fig. 3. — Uncalibrated color-magnitude diagram from //FUN data of a 6' square around 
OGLE-2003-BLG-238 that has been translated to put the centroid of the clump {triangle) 
at its known position [{V — /)o, -^o] clump = (1-00, 14.32). The unmagnified source (circle) is 
about 0.02 mag brighter and 0.22 mag redder than the clump. 



